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The motion of molecules across channels is critically important for understanding 
mechanisms of cellular processes. Here we investigate the mechanism of interactions 
in the molecular transport by analyzing exactly solvable discrete stochastic models. 
It is shown that the strength and spatial distribution of molecule/channel interac- 
tions can strongly modify the particle current. Our analysis indicates that the most 
optimal transport is achieved when the binding sites are near the entrance or exit of 
the pore. In addition, the role of intermolecular interactions is studied, and it is ar- 
gued that an increase in flux can be observed for some optimal interaction strength. 
The mechanism of these phenomena is discussed. 
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The role of biological channels is to support cellular processes by regulating fluxes of 
molecules and ions It is known that membrane protein pores are fast, efficient, selec- 

tive, and their functioning is robust with respect to strong non-equilibrium fluctuations in 
the cellular environment. These observations are especially surprising since in many cases 
the molecular translocation does not involve the use of metabolic energy or conformational 
changes 3(]. Understanding molecular mechanisms of transport is one of the most funda- 
mental problems in biological systems. 

When molecules enter into the channel their motion is slowed down mostly due to en- 
tropic barriers, and additional forces are needed to overcome these barriers. There are 
increasing experimental evidences that high efficiency, speed and selectivity of many bio- 
logical and artificial channels is a result of complex processes that involve molecule/pore 
and intermolecular interactions js-Sj. Recent high- resolution experiments on polypeptide 
translocations through protein nanopores 5|] started to probe the effect of molecule/pore 
interactions at the single-molecule level. It was found that changing the location of the 
binding site in the pore significantly modified the polypeptide fiux across the channel. How- 
ever, a comprehensive description of the role of interactions in the transport through the 
pores is still not available due to the strong biochemical and biophysical complexity of the 
translocation machinery and due to the lack of structural information 3|, |5|. 

To uncover mechanisms of molecular transport across the nanopores several theoretical 
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methods have been proposed [6H9|. The continuum models of the channel transport view 
the translocation as one- dimensional motion in an effective potential created by interactions 



with nanopores and with other molecules [7"]. A different approach utilized discrete-state 
stochastic models in which the translocation dynamics is analyzed as hopping between dis- 
crete binding sites in the channel js,!^!. Theoretical calculations show that both continuum 
and discrete approaches are closely related. Current theoretical models provide a reasonable 
description of some features of transport processes in the nanopores [7H9|. In this Letter, 
we analyze theoretically the effect of interactions on the molecular transport across chan- 
nels. Specifically, we address the question of how the translocation dynamics is modified 
by the strength and spatial distribution of the binding sites, and also by the intermolecular 
interactions. 

We consider transport of molecules in the nanopore as an effective one-dimensional motion 
along the discrete lattice of binding sites as illustrated in Fig. 1. There are N binding sites 
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in the channel, and concentrations of molecules to the left or right of the channel are equal 
to Ci and C2, respectively. The molecule can move into the channel from the left (right) with 
the rate uq = konCi {wq = konC2)] and the particle can move out of the channel with rates 
wi and un- see Fig. 1. In the nanopore the molecule at the site j (j = 1,2, ■• ■ , A^) can 
jump forward (backward) with the rate uj [wj). First, consider the situation when only one 
particle can be found in the channel. The probability to find the molecule at site j at time 
t is given by a function Pj{t), and the translocation dynamics is fully described by a set of 
master equations, 

dP(t) 

= Uj-iP,-i{t) + w;,+iP,+i(t) - {u, + w,)P,{t) (1) 

for j = 1, • ■ • , A^; while Po{t) = P/v+i(t) = 1 — Yl!^-i ^ji^) describes the completely empty 
channel at the time t Q]. We have shown earlier [81] that this model with binding sites 
can be solved exactly at t ^ 00 by mapping it into a single-particle random walk model 
on an infinite periodic lattice (with a period equal to + 1). Specifically, for the uniform 
channel with zero particle concentration to the right of the pore {wq = 0) the expression for 
the particle current is given by 

(A^ + 1) [u + yUo) 

where Uj = Wj = u {j = 1, ■ ■ ■ , N). To quantify the effect of interactions we assume that 
in one of the binding sites, say k, the particle interacts with the pore with potential e that 
differs from other sites. The case of £ > corresponds to attractive interactions, while 
negative e describe the repulsive binding site. The transition rates near the special binding 
site must satisfy the detailed balance conditions which lead to 

_ Uk-l u'f^ _ Uk 1 , . 

t — ^1 , — 1 l-jj 

where Mfc-i, Wfc, Wk and Wk+i correspond to the uniform channel without special interactions, 
and we define x = exTpie/knT). The corresponding explicit expressions for transition rates 



can now be written as 



= Uk-lX , Mfc = UkX , W'l, = WkX , W^+i = Wk+lX , (4) 

where the coefficient 6 (0 < 6 < 1) describes how the potential modifies the corresponding 
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transition rates [8|, 110| . Now fiux in the channel with the binding site at position k is equal 



4 



to 



Jk{N) 



u [2x-^ + iV - 1] + Mo 



UUq 

2{k - l)x-^ + x^-o + {N- k)x + ^^^^^ 



(5) 



The effect of interactions can be better understood by analyzing the ratio of particle currents, 

Jk{N) _ {N + l)[{u/uo) + f)] 



MN) {u/uo) [2x'' + AT - 1] + \2{k - l)x-' + x^-' + {N - k)x + ^^^^ -k+l 

' (6) 

The curves presented in Fig. 2 show how particle fluxes change depending on the position 
of the binding site. For attractive interactions the most optimal current is reached when 
the binding site is the last one {k = N), while it is better to have the repulsive site at the 
entrance {k = 1) to accelerate the transport. It can be shown rigorously from Eq. ([6]) that 
^'^^Q^^ > for positive e, and Jk{N) is always a decreasing function for negative e. These 
observations can be understood in the following way. Putting the attractive binding site 
near the exit increases the probability of finding the particle here, which leads to higher 
chances to complete the translocation by exiting to the right. The repulsive site at the 
entrance serves as a barrier for the particles that already passed it, lowering the probability 
of unsuccessful excursions without the translocation. These results are in agreement with 
single-molecule experiments on translocation of polypeptides js]. In these experiments the 
mutation in the biological nanopore that increased the molecule/pore interaction have led to 
faster transport when the mutation site was near the exit. It might also explain why so many 
biological channels have their binding sites at the entrance and at the exit positions since 



these distributions will optimize the overall fluxes Our results can be easily extended 
to more complex potential with several attractive and repulsive sites, and it can be shown 
that the most optimal flux is reached when repulsive sites cluster together near the entrance, 
while attractive sites tend to stay closer to the exit. 

The strength of interactions can also affect the flux through the nanopore as shown in 
Fig. 3, in agreement with previous theoretical findings for the channel-facilitated molecular 
transport [7|, l8| . For any set of parameters there is an optimal interaction strength e* that 
can be obtained from = 0, yielding 
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— + k-l 

Uo 



{1 - 9)x + {N - k)x 



i+e 



(7) 
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Specifically, for the most optimal site k = N (for attractive interactions) we have the 
following expression for the most optimal interaction strength, 



ksT In 
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— + N -I 



From Eq. ([7]) it can be shown that for 6* = we have e* 
binding site, while for 6 = 1 one can obtain 



1. 



-fcsTln 



2{^ + k-l 
N-k 



(8) 

-oo for any position of the 



(9) 



During the translocation across the channels more than one molecule can be found inside 
ihe nanopores and interactions between them might become important for the transport 
9|. Previous theoretical treatments considered the effect of the particle crowding [9], but 
assuming only hard-core exclusion interactions and neglecting correlations. To investigate 
explicitly the effect of intermolecular interactions we consider a specific N = 2 model without 
molecule/pore interactions, as specified above, but allowing more than one particle to be 
found in the pore. There is an energy cost associated with finding two particles next to each 
other. The configuration with 2 particles has an energy e, with e > {e < 0) describing 
attractive (repulsive) interactions. There are four possible configurations in the channel as 
plotted in Fig. 4. We label them as (i, j) with i,j = = 1) for the empty (occupied) 
site. It should be noted that the rate to enter the half-filled configuration ui and the exit 
rate from the fully occupied state U2 are related via the detailed balance, 

^ = ^x, (10) 

with X = exp{e/kBT). The case e = corresponds to the situation analyzed in Refs. |9|. 
This allows us to write explicit expressions, 

1 



Ui = UqX 



U2 = UX 



(11) 



where the coefficient < < 1 again specifies how the inter-particle interaction modifies 
these entrance and exit rates. We can define P{i,j;t) as the probability to find channel in 
the state at time t, and temporal evolution of the system dynamics can be found by 
analyzing corresponding master equations. Solving these equations at large times yields the 
expression for the molecular flux. 



J. 



uuo{u + ^x^ 



3u2 + ^(x + x'') + uuo(3 + 4) 



(12) 
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In the limit of e — oo only single molecules can be found in the channel and Eq. (fT2|) . 
as expected, reduces to Eq. (|2]) for = 2, Ji = 3^""°^^ - The effect of intermolecular 
interactions on the channel fluxes, as shown in Fig. 5, is rather complex. For = the flux 
is always a decreasing function of the interaction, and the single-particle transport is the 
most optimal. For 6 = 1 the trend is reversed: the stronger the interaction, the larger the 
molecular flux. However, for intermediate values of < ^ < 1 a non-monotonous dependence 
is observed with the flux reaching a maximum at some optimal interaction strength. The 
optimal interaction could be positive or negative depending on the parameters of the system. 
For attractive interactions the presence of the particle in the channel stimulates the entrance 
of another particle into the pore, but it slows down exiting of both particles from the channel. 
For repulsive interactions partially-filled channels serve as a barrier for the particle to enter, 
but simultaneously the entering particle accelerates the exiting of the particle inside. The 
combination of these processes explains the complex behavior in the channel with multiple 
particles. 

To summarize, we have investigated the effect of interactions on the molecular transport 
across channels. Using exactly-solvable discrete stochastic models, we have shown that the 
strength of the interaction, as well as the spatial distribution, are important parameters that 
can effectively control molecular translocations through nanopores. It was calculated that 
the largest particle current can be achieved when attractive sites are near the exit and/or 
repulsive sites are near the entrance. We have argued that the mechanism of how the 
interaction affect the transport across the channel is based on controling local concentration 
of particles in the channel. Attractive sites increase the probability to find the particles at 
these binding sites, while the repulsive sites work as barriers preventing particles already 
in the channel from moving back. Our theoretical picture agrees well with single-molecule 
experiments on translocation of polypeptides [s], and it might also explain distribution 
of binding sites in real biological channels In addition, we have studied the role of 

intermolecular interactions in the transport through nanopores. It was found that at some 
interaction strength the particle flux can be increased to reach the maximum level. The 
complex behavior could be explained by the fact that particles already in the channel might 
catalyze or inhibit the entrance into the channel of other particles. The presented theoretical 
model presents a theoretical framework for investigating complex transport phenomena in 
biological and artificial channels, and it might serve as a first step for further studies that 
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must include more realistic structural and biochemical information. 
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Figure Captions: 

Fig. 1. A schematic picture of the discrete stochastic model with N binding sites for the 
translocation through the pore. 

Fig. 2. (color online) The ratio of particle currents for different positions of the special 
binding site for the channel with = 10 binding sites. Circles are for £//cbT — 5, u/uq — 0.1 
and 9 — 0.5. Squares are for s/kBT — —5, u/uq — 0.1 and 9 — 0.5. Triangles are for 
s/kBT — 5, u/uq — 10 and 9 — 0.5. Diamonds are for s/kBT — 5, u/uq — 0.1 and 9 — 0.0. 

Fig. 3. (color online) The ratio of particle currents as a function of the interaction strength 
for the channel with = 10 binding sites. For the solid curve the parameters are u/uq — Q.1, 
k = 1 and 9 = 0.5. For the dotted curve the parameters are w/mq = 0.1, = 10 and 9 = 0.5. 

For the dashed curve the parameters are -u/mq = 0.1, /c = 10 and 9 = 0.8. For the dash-dotted 
curve the parameters are u/uo = 0.1, k = 5 and 9 = 0.9. 

Fig. 4. A general schematic picture for the channel with N = 2 binding sites and with 
intermolecular interactions. Open circles describe empty sites, while filled circles denote the 
occupied sites. 

Fig. 5. (color online) Ratio of the particle currents as a function of the intermolecular 
interaction for the channel with N — 2 binding sites. For the solid curve the parameters are 
u/uq — 0.1 and ^ = 0. For the dotted curve the parameters are u/uq — 0.1 and 9 — 0.5. For 
the dashed curve the parameters are u/uo — 0.1 and 9 — 1. For the dash-dotted curve the 

parameters are u/uq = 10 and 9 = 0.5. 
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Figure 1. A. Kolomeisky and K. Uppulury, Physical Review Letters. 
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Figure 2. A. Kolomeisky and K. Uppulury, Physical Review Letters. 
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Figure 3. A. Kolomeisky and K. Uppulury, Physical Review Letters. 




Figure 4. A. Kolomeisky and K. Uppulury, Physical Review Letters. 



